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The mechanisms producing strong coupling between electric and magnetic order in multiferroics are not 
always well understood, since their microscopic origins can be quite different. Hence, gaining a deeper 
understanding of magnetoelectric coupling in these materials is the key to their rational design. Here, we use 
ultrafast optical spectroscopy to show that the influence of magnetic ordering on quantum charge 
fluctuations via the double- exchange mechanism can govern the interplay between electric polarization and 
magnetism in the charge-ordered multiferroic LuFe 2 0 4 . 

Ferroelectric and ferromagnetic materials possess spontaneous electric and magnetic order, respectively, 
which can be switched by applied electric and magnetic fields. Multiferroics combine these properties in 
a single material, providing an avenue for controlling electric polarization with a magnetic field and 
magnetism with an electric field. These materials have been intensively studied in recent years, both for their 
fundamental scientific interest as well as their potential applications in a broad range of magnetoelectric 
devices 1 " 4 . 

Recently, the iron-based multiferroic LuFe 2 0 4 has attracted much attention because it exhibits magnetoelectric 
coupling close to room temperature 5 " 21 . The unique layered structure of LuFe 2 0 4 consists of double layers of Fe 
ions connected in a triangular lattice in the afr-plane (Figure 1(a)) 14 . The average valence of Fe ions is Fe 2 5+ , with 
Fe 2+ and Fe 3+ ions occupying equivalent sites in different layers with equal densities. The corresponding spin 
values are S = 2 (Fe 2+ ) and S = 5/2 (Fe 3+ ), with the spin structure shown in Figure 1(b) 1516 . A simple description 
based on nearest-neighbour interactions between Fe ions leads to the characterization of this material as a spin 
and charge frustrated system 517 . Bulk ferroelectricity was observed below the charge ordering temperature, T co ~ 
320 K, resulting in a spontaneous electric polarization that further increased upon the appearance of ferrimag- 
netic spin order below the Neel temperature, T N ~ 240 K 5 . Ferroelectricity in each bilayer is thus induced by 
electronic charge ordering, although the stacking of adjacent bilayers (i.e., in an antiferroelectric or ferroelectric 
arrangement) is still controversial 7,8,11,12,15,18 . Regardless, in each bilayer, the electric polarization P is coupled to 
the magnetic degree of freedom in LuFe 2 0 4 , but a comprehensive understanding of the mechanism underlying 
this magnetoelectric coupling has eluded researchers to date. Knowledge of this mechanism could potentially 
allow researchers to optimize both the strength of this coupling and its operating temperature to address the 
general goal of developing multiferroic materials with strong magnetoelectric coupling at room temperature. 

Theoretical studies have linked the magnetoelectric coupling in LuFe 2 0 4 to both thermal 19 and quantum 20,21 
charge fluctuations. In general, magnetic ordering can modify the effective hopping amplitude between two ions 
via the well-known double-exchange mechanism 20 " 23 , in which hopping is governed by the angle between the two 
core spins, as shown in Fig. 1(b). This change in the hopping amplitude will necessarily affect the quantum 
fluctuations of any charge ordered state of electronic origin. The corresponding change in charge ordering will 
necessarily modify P. In other words, if electronic charge ordering leads to a net electric polarization, the value of P 
should be modified by the presence of magnetic ordering (Fig. 1(b)). Moreover, if the magnetic ordering reduces 
(on average) the effective hopping amplitude, the corresponding suppression of quantum charge fluctuations 
leads to an increase of P below T N . Since the same mechanism should affect optically induced charge fluctuations, 
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Figure 1 | Charge and spin ordering in LuFe 2 0 4 above and below T N . (a) For T > T N , charge ordering results in a finite polarization P. The top layer is 
displaced from the bottom layer by an angle shown by the black straight dashed line 3 ; this is also shown clearly in Figure 3. Quantum fluctuations between 
Fe 2+ and Fe 3+ ions (depicted by blue dashed lines) can reduce Pby delocalizing charges, with an effective matrix element for hopping between two sites 
given by f 12 = til. Red arrows show the possible charge transfer routes between Fe 2+ and Fe 3+ , as defined in the text, (b) For T < T Ny ferrimagnetic spin 
ordering increases Pby decreasing the average hopping through the double exchange mechanism; t u = f when spins at both sites are ferromagnetically 
aligned (blue dashed lines), and t 12 = 0 when both spins are antiferromagnetically aligned (red dashed lines). 



we use femtosecond optical pump-probe spectroscopy, which has 
been extensively used to shed light on the properties of correlated 
electron materials 24 " 31 , to directly photoexcite and probe the 
Fe 2+ ^Fe 3+ charge transfer channel in LuFe 2 0 4 . Then, by varying 
the temperature T above and below T N , we can shed light on the role 
of these fluctuations in governing the coupling between spin and 
charge order in a single LuFe 2 0 4 bilayer, regardless of whether the 
bilayers are stacked ferroelectrically or antiferroelectrically. We find 
that the interlayer hopping matrix element describing these fluctua- 
tions depends strongly on their local core spin alignment via the 
double -exchange mechanism, making charge derealization (in real 
space, as shown in Fig. 1(b)) and hence the electric polarization 
extremely sensitive to the spin structure evolution over a broad 



temperature range. Therefore, although magnetoelectric coupling 
in various multiferroic materials has been studied using many dif- 
ferent techniques 32 , to the best of our knowledge, this is the first 
experimental evidence of magnetoelectric coupling mediated by 
the double- exchange mechanism in an insulator. 

Results 

We begin by developing a model for electronic hopping between two 
atomic sites, governed by the double exchange mechanism, which 
shows that the transition rate (i.e., charge transfer rate) between the 
ground and the excited state induced by an external driving electro- 
magnetic field is proportional to the effective hopping matrix, ty 2 , 
according to the Fermi Golden rule 33 , as described in Methods. 
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Figure 2 | Angle- dependent reflectivity measurements in LuFe 2 0 4 . (a) Imaginary part of the refractive index for p polarized excitation at incident angles 
of 25 and 45 degrees, (b) The data from (a), normalized to its value at 1.1 eV in order to compare the relative peak intensities at 1.5 eV. 
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Furthermore, the amount of delocalized charge 3q is proportional to 

Uj 2 , dq = (^) 2 ' where A is the energy difference between the ground 

and excited states. This simple observation establishes our ability to 
probe quantum charge fluctuations between two atomic sites using 
optical spectroscopy (see Methods for more detail). In LuFe 2 0 4 , 
these fluctuations are due to charge transfer between Fe 2+ and Fe 3+ 
ions (as revealed by optical spectroscopy 34 and band structure calcu- 
lations 19 ). We consider four different charge transfer channels in the 
bilayer crystal structure of LuFe 2 0 4 : interlayer charge transfer from 
the Fe 2+ rich bottom layer to the Fe 3+ rich top layer (£j) or from the 
top to the bottom layer (£|), and intralayer charge transfer within the 
top layer (E t _+) and within the bottom layer (E b ^), as shown in 
Fig. 1(a). We can gain insight on the relative energies of these dif- 
ferent charge transfer channels by considering the Coulomb energy 
between Fe ions in the Hamiltonian, 



where the pseudospin operators and Q z - are 1/2 or — 1/2 for Fe 3+ 
or Fe 2+ , respectively, and £ 0 > £ r and Ty are the permittivity of free 
space, the relative permittivity and the distance between sites i and 
j, respectively. Considering only the largest three interaction terms, 
we find that has the lowest excitation energy, E t _+ and E b ^ have 
intermediate excitation energies, and E± has the highest excitation 
energy. The interlayer transitions can be distinguished by noting that 
the bottom layer is rich in Fe 2+ while the top layer is rich in Fe 3+ . 
Therefore, it is clear that if the top layer has a positive charge density 
a > 0 per unit area, the bottom layer must have the opposite charge 
density, — a, to ensure charge neutrality. E± increases with a while E^ 
decreases, so it is reasonable to assume that E± > £| in LuFe 2 0 4 . We 
can also distinguish the intralayer transitions by noting that the 
configuration of in -plane oxygen ions around Fe ions in both layers 
leads to a higher in -plane charge transfer excitation energy for 
than that of E t ^. The optical conductivity measurements described 
in ref. 34 show two distinct charge transfer excitation channels at 
— 1.1 eV and —1.5 eV, which should thus correspond to £ T and E t _^, 
respectively. To further confirm this, we performed angle- dependent 
reflectivity measurements with a p polarized laser beam at different 
incident angles (25 and 45 degrees), which allowed us to extract the 
imaginary part of the refractive index (Figure 2(a)). Our data clearly 
shows two peaks at — 1 . 1 and 1.5 eV, in good agreement with the data 
from ref. 34, which is also plotted for comparison. It is worth noting 
that our experiments and those of ref. 34 were performed in different 
setups on samples grown by different groups; the close agreement 
between the data sets thus gives additional confidence in our results. 

Figure 2 supports our assignment of the two observed peaks to the 
energies E^ (1.1 eV) and E t _+ (1.5 eV), due to inter- and intra-layer 
charge transfer, respectively, as follows: As we increase the incident 
angle in our reflectivity measurements, we expect the absorption at 
£j to increase and E t _+ to decrease for p polarized light. This can be 
seen more clearly in Fig. 2(b), in which we normalized the data of 
Fig. 2(a) to the peak at 1.1 eV to facilitate the comparison of both 
peaks. This shows that the relative strength of the 1.5 eV peak, as 
compared to the 1.1 eV interlayer charge transfer peak, becomes 
weaker as we increase the incident angle, supporting our assignment 
of these two peaks. It is worth noting that both in our measurements 
and in the data of ref. 34, no spectral signatures corresponding to E b _+ 
and Ei were observed. This is likely because there are many different 
possible transitions that overlap at higher energies, which obscure 
the peaks corresponding to E b ^ and E±. Therefore, we used photon 
energies of 1.1 (£|) and 1.5 (£ f _>) eV in our experiments to examine 
inter- and intralayer quantum charge fluctuations in LuFe 2 0 4 . 

We propose that magnetic order and charge fluctuations in 
LuFe 2 0 4 are linked through the double exchange mechanism 20,21 , 
which leads to an effective hopping matrix element ty (see 



Methods) between the ions i and j that is determined by the angle 
Ojj between the spins S,- and S ; - : ty = tcos(6y/2) 23 . Within an Ising spin 
model, the hopping matrix element ty will be 0 (if the Fe 2+ and Fe 3+ 
core spins are antiparallel) or t (if they are parallel). Thermal fluctua- 
tions prevent any preferred spin orientation in the paramagnetic 
phase (T > T N ), so ty will be given by its average value of t/2 
(Fig. 1(a)) (similarly, ty 2 = 0 or t 2 , respectively, and thus its average 
value will be f/2 for this transition). However, in the magnetically 
ordered state (T < T N ), ty = f if the nearest neighbour spins are 
aligned ferromagnetically (FM), while ty = 0 if they are antiferro- 
magnetically (AFM) aligned (Fig. 1(b) and Figure 3). This can then 
be applied to the four charge transfer channels discussed above (con- 
sidering electron hopping from a given Fe 2+ ion to its nearest neigh- 
bour Fe 3+ ions) to find the total steady- state absorption (propor- 
tional to ty 2 (see Methods)) for each transition above and below 
T N . We find that the development of magnetic order does not affect 
the total ty 2 for E^ and E t ^. This can straightforwardly be seen by 
considering the E^ transition (Fig. 1(a)), where the spins do not have 
a preferred direction above T N . Because there are two possible charge 
transfer channels for E^ (Fig. 1(a)), we obtain a total transition ampli- 
tude of t 2 /2 X 2 = f above T N . Below T N , there is only one possible 
transition with an amplitude of t 2 (Fig. 1(b)), so the transition ampli- 
tudes above and below T N are the same. In the same manner, if we 
consider the E t _+ transition above T N (Fig. 1(a)), there are six charge 
transfer channels, giving a transition amplitude of t 2 /2 X 6 = 3t 2 . 
Below T N , there are only three possible transitions, each with ampli- 
tude t 2 , so the transition amplitudes above and below T N are the 
same, as for the interlayer transition. 

This argument thus verifies that there is no change in the steady- 
state absorption for both transitions examined here (E t _+ and E^) 
between T > T N and T < T N . This is also clear from the data in 
ref. 34, which shows that the absorption for the E^ and E t _+ transi- 
tions does not significantly change with temperature. In contrast, the 
total ty 2 is reduced through the double exchange mechanism for E b ^ 
and £j below T N ; therefore, as described above in the introduction, 
the corresponding suppression of quantum charge fluctuations 
should increase P. A simple test of the influence of the double 
exchange mechanism on magnetoelectric coupling in LuFe 2 0 4 
would thus be to track the expected changes in the steady-state 
optical absorption at E b ^ and E± as the temperature Tis varied across 
T N . However, the conductivity data in ref. 34 does not show any 
distinct spectral signature that could be attributed to these charge 
transfer energies, likely due to the presence of many overlapping 
transitions in this energy range as described earlier. 

One can circumvent this limitation by performing a non-equilib- 
rium experiment; i.e., photoexciting the known E^ and E t _^ transi- 
tions and examining the resulting changes in the optical absorption 
at the £ T transition (which is proportional to P). Photoexcitation 
changes the charge configuration, which in turn changes both the 
total ty 2 and the £j charge transfer energy through equation (1), 
leading to a transient change in the reflectivity that can be measured 
in our experiments (see Methods for more detail). The development 
of magnetic order below T N can then further modify the absorption 
at £ j through the double exchange mechanism in the same manner as 
in the steady state, potentially causing an additional change in the 
photoinduced reflectivity. In short, our ultrafast optical experiments 
(described in more detail in Methods) allow us to examine the effect 
of photoexciting either intralayer or interlayer charge transfer transi- 
tions (in effect externally driving charge fluctuations) on the inter- 
layer charge transfer energy as a function of temperature. Tuning the 
temperature above and below T N then allows us to determine the 
effect of magnetic ordering on charge fluctuations. In this way, we 
can test if these fluctuations are indeed responsible for the magneto - 
electric coupling measured in LuFe 2 0 4 . 

Figure 4(a) shows the temporal profile of the normalized photo - 
induced reflectivity change, ARIR(t), in LuFe 2 0 4 at several different 



SCIENTIFIC REPORTS | 3 : 2654 | DOI: 1 0.1 038/srep02654 



3 



(a) (b) 



After l.leV excitation After 1.5eV excitation 




Figure 3 | Possible interlayer and intralayer charge transfer transitions for T< T N . Both figures show the spin structure of LuFe 2 0 4 for T< T N . Big and 
small solid circles correspond to atoms in the top and bottom layers, respectively. The black arrows show the local spins along the c axis at 220 K, from refs. 
15 and 16. Red solid and broken arrows represent allowed and forbidden interlayer (a) and intralayer (b) charge transfer channels, respectively, when 
considering the double exchange mechanism. After 1.1 eV excitation, one electron from an Fe 2+ ion in the bottom layer moves to an Fe 3+ ion in the top 
layer, while after 1.5 eV excitation, one electron from an Fe 2+ site in the top layer moves to an Fe 3+ site in the same layer. 



temperatures for a degenerate 1.1 eV pump-probe measurement. 
Immediately after photoexcitation, ARIR decreases to its minimum 
value within <0.5 picoseconds (ps) (Figure 4(b)) and returns to 
equilibrium while exhibiting coherent acoustic phonon oscillations 
with a period of —40 ps (Fig. 4(a)). These coherent phonons are 
generated by the dynamic stress on the sample induced by absorption 
of the pump pulses and are commonly observed in ultrafast optical 
experiments on correlated electron materials 30 , as well as non-corre- 
lated materials described by a simple band structure 35 . In this paper, 
we will focus on the variation of the maximum amplitude of the 
transient reflectivity signal (AR/R max ) with temperature, which gives 
insight into the influence of intra- and interlayer charge fluctuations 
on magnetoelectric coupling in LuFe 2 0 4 . 



Figure 4(c) depicts ARIR max (1.1 eV) after 1.1 eV photoexcitation 
as a function of temperature near T N . As described above, photoex- 
citation changes the charge distribution in LuFe 2 0 4 , altering the 
energy required to transfer a charge from the bottom to the top layer 
(£|) (and thus the absorption/reflectivity probed at 1.1 eV) (Figure 3 
(a)). We calculated AEp the difference in the interlayer charge trans- 
fer energy before and after photoexcitation, for T < T N and T > T N , 
using equation (1) (and using s r = 2 from ref. 34) (see Methods). This 
calculation reveals that there is no change in the ARIR max signal as T 
is varied across T N , which is consistent with our experimental obser- 
vation (Fig. 4(c)); fundamentally, although the specific allowed 
charge transfer transitions after photoexcitation do change across 
T N due to the double exchange mechanism (Fig. 3(a)), the net effect 
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Figure 4 | Temperature-dependent transient reflectivity change after 1.1 eV photoexcitation. (a) Temperature-dependent transient reflectivity 
change after photoexciting and probing the Fe 2+ — > Fe 3+ interlayer charge transfer excitation at 1.1 eV. The open circles are the experimentally obtained 
data points and the solid lines are the results of fitting the data with exponential and oscillating terms, (b) The data from (a) at early times, (c) The 
amplitude of the negative peak as a function of temperature. The red vertical line shows the magnetic transition temperature, T N ~ 240 K, and the blue 
horizontal line is parallel to the x axis, which indicates that there is no significant variation in AR/R max (1.1 eV) across T N . We note that the high signal- 
to-noise ratio of our experiment resulted in error bars that are comparable to the size of the data points and are therefore not shown. 
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Figure 5 | Temperature-dependent transient reflectivity change after 1.5 eV photoexcitation. (a) Transient reflectivity change at 1.1 eV for 
temperatures of 200 and 300 K after photoexciting the Fe 2+ — » Fe 3+ intralayer charge transfer excitation at 1.5 eV. (b) Peak amplitude of the transient 
reflectivity change near T N after 1.5 eV excitation. 



of magnetic ordering on the resulting absorption change is very small 
and thus undetectable in our experiment. The negative sign of the 
signal is also expected since photoexcitation reduces the absorption 
at 1.1 eV. 

Figure 5(a) shows the transient reflectivity change for the 1.1 eV 
interlayer transition after photoexciting the intralayer (E t _+) charge 
transfer channel at 1.5 eV. The time- dependent dynamics are similar 
to those observed after photoexciting Ep but the variation of ARI 
Rmax wim temperature is quite different (Figure 5(b)); in particular, a 
significant increase in the amplitude is clearly observed as the tem- 
perature rises above T N . As described above, the steady state absorp- 
tion for the 1.1 eV interlayer transition does not change across 
T N ; however, the probe absorption at this transition can change after 
1.5 eV photoexcitation as Tis varied across T N . To understand this, 
we calculated the effect of the photoinduced intralayer charge trans- 
fer at 1.5 eV on in the same manner as described above for £ j (see 
Methods). Our calculation shows that the ratio of the maximum 
photoinduced change in reflectivity between T < T N and T > T N 
is —0.9, which agrees very well with our experimental results (—0.87) 
(Fig. 5(b)). 

Discussion 

From this experimental observation, we deduce that the ferrimag- 
netic order influences the fluctuations of the charge ordered state that 
is responsible for the electric polarization in LuFe 2 0 4 through the 
double exchange interaction. This result indicates that the interplay 
between charge fluctuations and magnetic ordering can result in 
magnetoelectric coupling at the Neel temperature. Finally, it is worth 
noting that this mechanism will generate an electric polarization in 
each bilayer, regardless of whether the ground state consists of layers 
stacked with ferroelectric or antiferroelectric order. 

In summary, we used femtosecond optical pump-probe spectro- 
scopy to investigate the role of the double exchange mechanism in 
the magnetoelectric coupling observed in LuFe 2 0 4 . Our experi- 
ments revealed that optically induced charge fluctuations are affec- 
ted by magnetic order in a manner that is consistent with this 
mechanism. Importantly, this result opens an alternative route for 
finding strong magnetoelectric effects: charge ordering in transition 
metal oxides can naturally lead to electric polarization that is 
coupled to the magnetic degree of freedom via the double-exchange 
interaction. 



Methods 

Experimental setup and sample preparation. Our femtosecond optical pump-probe 
spectroscopy system is based on a 75 femtosecond (fs), 250 kHz repetition rate 
Ti:sapphire amplified laser system operating at 800 nm (—1.5 eV) and seeding an 
optical parametric amplifier (OP A) that allows us to tune the photon energy. A delay 
line allows us to vary the optical path difference between the pump (200 urn 
diameter) and probe (100 urn diameter) beams, which are then focused to the same 
spot on the sample. Temperature-dependent transient reflectivity changes were 
obtained in reflection with cross-polarized pump and probe beams in a > 10 : 1 power 
ratio (the pump fluence is 76 uj/cm 2 , which photoexcites —0.007 electrons/unit cell), 
incident at an angle of less than 10° to the hexagonal c axis of the crystal. The probe 
photon energy was 1.1 eV in all experiments (examining E^), and the pump photon 
energies were 1.1 and 1.5 eV (photoexciting and E t _^, respectively). We emphasize 
that the interlayer Fe-Fe transition dipole is not oriented perpendicular to the surface 
(and thus not aligned with the c-axis since Fe ions in the top layer are not located 
above Fe ions in the bottom layer, but instead above O ions) (Figs. 1 and 3) 3 ' 6 ' 13 ' 15 ' 16 , so 
our near normal incident 1 . 1 eV pulses can sensitively photoexcite and probe changes 
in its absorption. It is also worth mentioning that pump fluence-dependent 
measurements revealed that only the amplitude of the ARIR signal changes linearly 
with fluence, with no changes in the measured dynamics, for both 1.1 and 1.5 eV 
excitation. Furthermore, at the fluence used here, the maximum temperature increase 
due to laser heating is <8 K, which should not significantly affect the measured 
dynamics, and the sample completely recovers in the 4 us time interval between 
amplifier pulses. Finally, the LuFe 2 0 4 single crystal used in this study was grown by 
the floating zone method as described in ref. 15, with its surface normal to the c-axis. 

Theoretical background. Quantum charge fluctuations originate from hopping of an 
electron between two spatially separated potential minima. When an electron is 
localized in one potential well, the system becomes electrically polarized 20 ' 21 . If we 
consider two energetically non-degenerate states localized at two different lattice sites 
and include the double exchange interaction in the system, the Hamiltonian can be 
expressed as 

H° = e 1 C+ q + s 2 C+ C 2 - t 12 q + C 2 - 1 21 C + Ci , (Ml) 

where £; is the atomic energy, and Q are creation and annihilation operators at the 
z'th site (i = 1,2), and ty = is a effective hopping matrix element accounting for the 
hopping between sites. This matrix element is governed by the double exchange 
mechanism, which relates the probability of an electron hopping between two atoms 
to the angle 9^ between the local core spins S f and S ; 22 ' 23 . This mechanism has most 
frequently been used to explain the metallic properties of colossal magnetoresistive 
manganites 22 , but also applies here to LuFe 2 0 4 , since the ferrimagnetic spin order 
existing below T N influences electron hopping (and therefore charge fluctuations of 
the charge-ordered state), which can in turn change the dielectric properties of the 
system. In other words, since the electronic charge is never completely localized in 
insulators, the degree of derealization depends on the effective hopping amplitude 
given by the double exchange mechanism. Since magnetic ordering suppresses this 
hopping amplitude for the E b ^ and E± transitions, we expect electrons in LuFe 2 0 4 to 
be more localized, stabilizing charge order. 

For small ty values (ty <C A, where s 2 — z x = A and t 12 = t 2 i) and e 2 > most of 

the charge will be localized at site 1, with a small fraction dq = ( ^ J of delocalized 
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charge remaining at site 2. Because the electric polarization P is proportional to the 
difference of charge densities between sites, p 2 ~ Pi (where p 2 and p x are the electron 
densities at sites 2 and 1, respectively), any change in the delocalized charge at site 2 
causes a change in P. 

When an external electromagnetic field (£ 0 coscof) is applied to the system, it will 
introduce a small perturbation H 1 = exE 0 coscdt into the Hamiltonian, inducing a site- 
to -site transition (where e is the electron charge and x is the distance between the two 
sites). We can use Fermi's golden rule 33 to calculate the probability of transitions 
between both sites (corresponding to quantum charge fluctuations), which is 
found to be proportional to the delocalized charge on site 2, Sq, through 

H 2 l2 = \( x ¥ 2 \H 1 \ x ¥ l )\ 2 -C-jf = 5q (where 1^) and |¥ 2 ) are the ground and 

excited states of H°, respectively). Note that the transition rate is proportional to the 
extent of charge derealization in the ground state of the system. Since the photo - 
induced change in reflectivity at the absorption peak is proportional to changes in the 
absorption under the conditions of our experiment 31 , which, in turn, is proportional 
to H\ 2 and 8q, we can relate our transient reflectivity measurements to the amount of 
delocalized charge and thus the polarization P. This then establishes that we can use 
our ultrafast optical measurements to reliably photoexcite and probe quantum charge 
fluctuations in LuFe 2 0 4 . 

We calculated AEp i.e. the pump-induced difference in the interlayer charge 
transfer energy, using equation (1) and only considering four ions: two excited by the 
pump and two examined by the probe after photoexcitation (Figure 3 in our manu- 
script). This was done by exchanging a Fe 2+ and a Fe 3+ ion either between the bottom 
and top layers (corresponding to absorption of a 1.1 eV pump photon) (Figure 3(a)) 
or between two sites in the top layer (corresponding to absorption of a 1.5 eV pump 
photon) (Figure 3(b)). We then obtained AE^ based on the new Coulomb energy after 
exchanging Fe 2+ and Fe 3+ ions. This was then used to calculate the change in 
absorption as described above, from which we calculated the variation of ARIR max 
with temperature after both 1 . 1 eV and 1 .5 eV photoexcitation for comparison to our 
experimental data. It is worth noting that if either of the ions that absorb a pump 
photon is involved in the subsequent absorption of a probe photon, we find that the 
resulting AE^ is much larger than the probe bandwidth (—13 meV) and thus does not 
contribute to the observed absorption change. Finally, including more than 47 Fe 2+ 
electrons in this calculation resulted in an insignificant reflectivity change, since there 
is almost no change in due to Fe 2+ ions far from the Fe 2+ and Fe 3+ sites that 
participate in the photoinduced transition. 
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